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ABSTRACT 

A  method  based  on  rational  mechanics  is  presented  for  predicting 
hydrodynamic  forces  and  moments  on  a  hull  of  arbitrary  geometry  as  a  re¬ 
action  to  the  propeller-induced  velocity  field.  This  method  assumes  that 
the  hull  surface  can  be  represented  by  a  distribution  of  sources.  The 
numerical  procedure  is  adapted  to  the  CDC  7600  or  Cyber  176  digital  com¬ 
puter  which  furnishes  the  strengths  of  the  source  distribution  over  the 
hull  and  the  vertical  component  of  the  propeller-induced  hull  force  at 
blade  frequency.  The  program  is  executed  for  two  widely  different  cases 
for  which  experimental  measurements  are  available:  1)  a  spheroidal  head 
in  the  presence  of  a  3-bladed  propeller  operating  in  uniform  inflow,  and 
2)  a  Series  60  hull  model  (V-form  stern)  driven  by  a  4-bladed  propeller. 
Theoretical  predictions  show  agreement  with  corresponding  experimental  measure¬ 
ments  within  11%  for  the  body  of  revolution  and  27%  for  the  surface  ship. 

A  disadvantage  of  this  approach  is  the  excessive  computing  time 
required  not  only  because  the  oscillatory  nature  of  the  propeller  velocity 
field  dictates  that  the  hull  be  divided  into  a  large  number  of  small 
quadrilaterals  but  also  because  the  representation  of  the  hull  surface  by 
sources  requires  the  evaluation  of  three  components  of  the  normal  velocity 
induced  by  the  propel  le-",  for  all  possible  combinations  of  frequencies  of 
the  loading  and  propagation  functions. 

A  new  method  has  lately  been  devised  which  replaces  the  sources  by 
a  distribution  of  doublets,  thus  replacing  the  three  components  of  normal 
velocity  by  the  negative  of  the  velocity  potential  which  has  no  directiv¬ 
ity.  Preliminary  results  show  a  substantial  reduction  in  computing  time 
and  improvement  in  the  accuracy  of  the  results  as  well.  The  calculated 
theoretical  predictions  for  the  same  quadrilateral  distribution  now  show 
an  agreement  with  the  measurements  within  2%  for  the  body  of  revolution 
and  10%  for  the  surface  ship. 

It  should  be  realized,  however,  that  Lewis's  measurements  involved 
the  combined  effect  of  a  force  and  a  couple,  in  contrast  to  the  theo¬ 
retical  prediction  which  gives  only  a  vertical  force.  Thus,  a  one-to-one 
comparison  is  currently  not  possible. 
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angular  veloci ty/f ree  stream  velocity 
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propeller  diameter 
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blade  thickness  distribution 
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i  ndex 
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&a 

A(R,(  ) 

\ 


velocity  component 

resulting  velocity 
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rectangular  coordinate  system  of  hull 

rectangular  coordinate  system  of  propeller 

cylindrical  coordinates  of  point  in  space 
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angular  chordwise  location  of  loading  point  on  propeller  blade 

result  of  chordwise  (0a)  integration  for  assumed  chordwise  mode 
(see  Reference  5) 

positive  integer 


§,Tj,£  rectangular  coordinates  of  source  distribution  on  hull 

§,p,0Q  cylindrical  coordinates  of  loading  point  on  propeller 

Pj-  mass  density  of  fluid 

source  density  on  hull 


a  angular  measure  of  skewness  of  propeller  blade 

t  variable  of  integration 

§  total  velocity  potential 

diffraction  potential  of  hull  boundary 
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cpD  diffraction  potential  of  hull 
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ai  angular  frequency 


VI  i  i 


» 


R-1979 


» 


» 


> 


* 


» 


r 


t 


v 


¥ 


INTRODUCTION 

One  of  the  primary  considerations  of  naval  architects  in  designing 
propeller-stern  configurations  is  minimization  of  propeller-induced  hull 
vibration.  One  characteristic  which  determines  the  selection  of  any 
particular  design  over  another  is  the  expected  level  of  vibratory  excita¬ 
tions  at  integer  multiples  of  blade  frequency. 

A  method  based  on  rational  mechanics  has  been  evolved  at  Davidson 
Laboratory  for  predicting  the  forces  and  moments  generated  on  the  shaft 
and  on  the  hull  surface  for  arbitrary  propeller  hull  geometries  and  their 
given  wakes.  It  couples  computat ional  methods  for  predicting  unsteady 
propeller  blade  loadings,  and  the  field  point  inductions  due  to  loadings 
and  blade  thickness  effects,  to  a  computational  method  for  the  hydrodynamic 
reactions  of  a  hull  of  specified  geometry. 

The  forces  generated  by  the  propeller  which  excite  ship  structural 
vibration  are  composed  of  those  generated  by  the  fanning  action  of  the 
blades  on  the  neighboring  hull  and  those  generated  on  the  shaft,  namely, 
the  bearing  forces  and  moments.  The  surface  forces  depend  on  all  shaft 
frequencies  of  the  blade  loading,  while  the  blade  frequency  bearing 
forces  and  moments  are  caused  by  various  loadings,  some  at  the  blade  fre¬ 
quency  of  interest  and  others,  one  above  and  one  below  the  blade  frequency. 
The  propeller  operates  in  a  spatially  varying  inflow  of  axial,  radial  and 
tangential  components  of  velocity.  Within  the  assumptions  of  linearized 
propeller  theory,  the  blades  respond  to  the  variations  in  the  axial  and 
tangential  velocities  only.  As  the  blades  traverse  the  wake,  they  respond 
to  the  wake  harmonics  in  such  a  manner  that  the  thrust  and  torque  are  due 
to  integer  multiples  of  the  blade  frequency  wake  harmonics,  and  the  bear¬ 
ing  forces  and  moments  are  due  to  frequencies  one  unit  above  and  one  unit 
below  the  integer  multiple  of  blade  f  requency  (iiN+1  and  £N-1  times  W , 
where  N  is  the  number  of  blades,  Q  is  the  shaft  angular  velocity,  and  £ 
is  the  multiple  of  interest.)  These  forces  are  the  sum  of  those  acting 
on  each  blade.  However,  each  blade  responds  to  all  the  harmonics  of  the 
circumferential  variation  of  the  wake.  Each  of  these  produce  a  contribu¬ 
tion  to  the  radiated  or  induced  field  of  the  propeller  at  integer  multiples 
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of  blade  frequency.  In  addition,  the  finite  thicknesses  of  the  blades 
simulated  by  source  distributions  produce  a  blade-frequency  radiation  which 
is  also  accounted  for  in  the  present  analysis.  In  the  method  presented 
here,  the  components  of  loading  due  to  the  wake  and  to  the  thickness  are 
considered  separately,  and  synthesized  at  the  end.  This  procedure  is  both 
convenient  and  appropriate  because  the  theory  is  linear. 

The  objective  of  the  present  investigation  is  to  develop  a  theoretical 
approach  to  evaluate  the  propeller-induced  vibratory  hull  forces.  Details 
of  the  development  of  the  mathematical  models  applied  to  the  unsteady 
propel ler/hul 1  interaction  problem  are  presented  in  this  report.  They  are 
followed  by  descriptions  of  the  numerical  analysis  procedure  and  the  coded 
algorithm.  The  vertical  hull  force  is  calculated  for  two  cases  for  which 
experimental  data  are  available.  The  agreement  obtained  for  these  first 
calculations  tends  to  substantiate  the  approach  taken  to  solve  the  unsteady 
propel ler/hul I  interaction  problem;  however,  it  is  recognized  that  addi¬ 
tional  checking  is  advisable  before  this  tool  is  to  be  applied  on  a  regular 
basis.  Comparisons  with  other  cases  available  in  the  literature  are  in 
progress. 

The  present  approach  requires  distribution  of  sources  over  the  hull 
surface,  the  strengths  of  which  are  determined  through  the  solution  of  a 
Fredholm  integral  equation  of  the  2nd  kind  in  which  the  known  function 
represents  the  propeller-induced  velocity  normal  to  the  hull  surface.  In 
an  investigation  currently  underway  at  Davidson  Laboratory,  the  sources  are 
replaced  by  doublets  whose  strength  are  determined  by  the  same  integral 
equation  but  with  the  known  function  replaced  by  the  negative  of  the  pro¬ 
peller  velocity  potential  function,  which  has  no  directivity  and,  hence, 
there  is  only  a  single  value  at  each  point  of  the  hull  instead  of  three 
components  of  the  normal  velocity.  This  should  reduce  computer  time  con- 
s iderably. 

This  research  was  sponsored  by  the  Maritime  Administration  under 
Contract  N00014-76-C-0862  and  administered  by  the  David  Taylor  Naval  Ship 
Research  and  Development  Center. 
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BACKGROUND 

In  1965,  Breslin  and  Eng'  were  the  first  to  attempt  to  evaluate  the 

prope 1 1 er- i nduced  hull  surface  forces  (i.e.,  the  fanning  forces)  by  com- 

2  3 

bining  an  extended  form  of  the  Hess-Smith  program  with  a  program  developed 
at  Stevens  Institute  of  Technology  to  calculate  the  propel  1 er- i nduced  veloc¬ 
ity  field.  The  latter  program  was  based  on  representing  the  propeller  blades 
by  lifting  lines  and  was  for  a  propeller  operating  in  uniform  inflow  having 
thus  steady  state  (time  invariant)  blade  loading  only.  A  numerical  procedure 
was  developed  for  the  IBM  7090  digital  computer  which  had  to  be  abandoned  as 
impractical  because  of  the  excessive  time  required  by  computers  of  that  era. 

4 

More  recently,  Vorus  developed  a  method  which  seeks  to  avoid  solving 
the  diffraction  or  scattering  problem  posed  here  by  generating  the  hull  in 
the  presence  of  the  nonuniform  onset  flow  from  the  propeller.  His  method, 
which  uses  a  reciprocity  theorem,  produces  a  formula  for  a  particular  hull 
force  component  in  terms  of  integrations  over  the  propeller  blades  of 
quantities  (approximated  by  asymptotic  expansions)  induced  at  the  blades  by 
unit  velocity  motion  of  the  hull  in  the  plane  of  the  desired  force.  Thus 
to  determine  three  forces  and  three  couples,  six  boundary-value  problems 
must  be  solved.  In  contrast,  the  method  advanced  here  requires  the  solution 
of  only  one  boundary-value  problem  and  all  forces  and  moments  are  obtained 
by  integrations.  Moreover,  by  using  the  Hess-Smith  program,  the  arbitrary 
hull  shape  is  exactly  taken  into  account  whereas  Vorus  uses  a  strip  theory 
to  model  the  hull;  strip  methods  tend  to  be  inaccurate  near  the  stern  of 
three-dimensional  bodies. 

Since  the  first  attempt  to  compute  the  unsteady  propel ler/hu 1 1  inter¬ 
action  problem  by  Breslin  and  Eng,'  systematic  investigations  have  been 
undertaken  to  improve  propeller  representations  and  to  treat  more  realistic 
operating  conditions.  An  unsteady  lifting  surface  theory  has  been  developed 


Superior  numbers  in  text  matter  refer  to  similarly  numbered  references 
listed  at  the  end  of  this  report. 
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and  coded  for  computing  the  steady  and  unsteady  propeller  blade  forces  and 
moments  taking  into  account  finite  chord  lengths,  helicoidal  geometry, 
thickness,  mutual  blade  interaction,  and  spatial  nonuniformity  of  the  in¬ 
flow  into  the  propeller  disc.  A  detailed  documentation  of  the  computer 
program  developed  on  a  CDC  6600  high-speed  digital  computer  has  been 
issued. The  program  furnishes  a  procedure  to  compute  the  mean  (time 
averaged)  and  vibratory  blade  pressure  distributions  together  with  the  in¬ 
tegrated  hydrodynamic  forces,  moments  and  blade  bending  moments  about  any 
face-pitch  line.  Furthermore,  it  provides  the  necessary  information  for 
the  study  of  inception  of  cavitation  as  well  as  for  the  study  of  blade 
stress  analysis  which  can  be  performed  by  combining  the  propeller  program 
with  the  STARDYNE-CDC  program  (or  the  like)  based  on  a  finite  element  model 
of  the  blade  structure.^ 

Based  on  the  unsteady  lifting  surface  theory  of  marine  propellers 

developed  at  the  Davidson  Laboratory,  theoretical  methods  and  computer  codes 

have  been  developed  to  compute  the  mean  and  vibratory  velocity  and  pressure 

fields  induced  by  a  noncav i ta t i ng  propeller  operating  in  both  uniform  and 

nonuniform  inflow  fields.  Both  loadinq  and  thickness  effects  have  been  taken 

8  9 

into  account  in  evaluating  the  propeller-induced  velocity  field.  ’ 

In  the  present  investigation  a  form  of  the  extended  Lagally  Theorem 
is  developed  and  programmed  to  compute  the  propeller-induced  hull  forces. 

The  velocity  field  program  based  on  the  unsteady  lifting  surface  theory  was 
exercised  extensively  before  it  was  merged  with  the  Lagally  Theorem  code  and 
the  Hess-Smith  code  to  provide  a  self-contained  computer  code  for  calculating 
the  propeller-induced  vibratory  forces.  The  details  of  this  code  are  des¬ 
cribed  in  the  sections  to  follow. 
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PROPELLER-INDUCED  VIBRATORY  FORCES 


In  this  section  a  derivation  of  the  extended  Lagally  Theorem  is  given 
which  will  be  utilized  for  calculating  the  hull  surface  forces.  The  vector 
equation  of  motion  of  an  incompressible  inviscid  fluid  is  taken  to  be: 


||  +  (u  •  v)u  =  -  i-  v  pT  (1) 

where  H  is  the  velocity,  p^  is  the  total  pressure  (including  perturbation), 
and  pj.  the  mass  density. 

The  velocity  field  satisi  js  the  continuity  equation 

V  •  “  =  0  (2) 

and  upon  joining  the  condition  of  i rrotat ional i ty ,  the  total  velocity  field 
may  be  written  in  terms  of  a  potential  defined  by 

u  =  vi  (3) 


From  Eq.(2)  it  is  observed  that  V2$  =  0.  Substituting  Eq.(3)  into  Eq.(l) 
and  integrating  once  yields 


M  1 

at  2 


v§  • 


W 


which  is  a  form  of  the  unsteady  Bernoulli  equation. 

The  total  velocity  potential  for  the  flow  under  consideration  may  be 
written  as 

$(a,t)  =  Ux  +  tps(><)  +  cpp(x,t)  +  9Q(x,  t) 
or 

§(x,t)  =  Ux  +  cps(x)  +  cp(x,t)  (5) 


where  Ux  is  the  potential  of  the  onset  flow  from  &  —  -<*>,  cp<.(x)  is  the 
diffraction  potential  of  the  hull  boundary  in  response  to  the  potential 
Ux,  cpp(x,t)  is  the  propeller  radiated  potential,  and  cpp , t)  is  the 
diffraction  potential  of  the  hull  ir.  response  to  the  onset  flow  described 
by  cpp(x,t)  which  is  imposed  by  the  propeller.  For  convenience,  the  unsteady 
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contributions  to  the  total  potential  are  combined  into  a  single  function 
cp(£,t).  Simple  harmonic  time  dependence  will  be  assumed,  viz., 

?(x,t)  =£  cM*)e,'eNU)t  (6) 

£= 0  x 

in  which 

v3^  =  0  (7) 

The  boundary  conditions  at  £  -*  -«  and  Z  — are  given  by  vcp^  —  0. 
q>(-°°,t)  is  taken  to  be  zero. 

The  boundary  condition  of  the  hull  surface  is  given  by 

n  •  Vc ?.(£)  =0  x  on  S  (see  Figure  1)  (8) 

where  n  is  taken  to  be  the  outward  unit  normal  (into  the  main  body  of  fluid). 
The  free  surface  condition  at  Y=0  is  taken  to  be 

9£(i)  =  0  Y=0  (9) 


On  substituting  Eq.(5)  into  Eq.(4),  defining  u$  =  Ui  +  V<+>s  (*)  [which 
is  the  total  time-independent  convective  velocity  around  the  hull],  and 
linearizing,  the  following  equation  of  motion  for  the  unsteady  perturbations 
is  obtained 


+  u  .  vcp  =  -  £- 
dt  ~s  T  pf 


(10) 


where  p  is  the  part  of  p^  due  to  cp  alone. 

The  velocity  potential  cp(=  cp^  +  cp^)  satisfies  Laplace’s  equation 

y3cp  =  V^cpp  =  V2cpo  =  0 

and  the  following  boundary  conditions 


9p  =  90  =  0 


5n 


=  0 


on  Y=0  (the  free  surface) 


on  S 


(11a) 

(lib) 


where  Y  and  S  are  defined  in  Figure  1. 

The  hull  surface  force  may  be  found  by  integrating  the  surface  pres¬ 
sures  (Eq.10)  over  the  hull  surface  and  its  negative  reflection  on  the 


t 

Ji 
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waterplane.  In  index  notation  the  components  of  hull  forces  due  to  the 
propeller  imposed  flow  field  are 


■ 


3(VV 

dt 


dX. 


ax. 


ST  ds"pf  y  Lv*j  SxT  VVVJ  57T  ds  (,2) 


where  Vs.  is  the  j-component  of  the  resulting  time-independent  velocity 
(as  defined  previously). 


In  the  present  investigation,  the  body  or  hull  surface  is  modeled  by 
a  distribution  of  sources  and  negative  source  images  reflected  about  the 
(X,Z)-plane  (to  model  the  free  surface  as  a  high  Froude  number  problem  —  see 
Figure  1).  Therefore,  the  diffraction  potential  may  be  written  as 

cpD(X,Y,Z;t)  =  -Jf  o(S,TU;t)  (£  -  i  )dS*  (13) 

S 1  1 

where  cr(§,Tl,G;t)  is  the  source  density, 

1 

R  =  |(X-§)2  +  (Y-T))2  +  (Z-G)2}2 
Ri  =  |(X-§)3  +  ( Y+7]) 2  +  (Z-C)2i2 

and 

S'  =  S 


The  temporal  part  of  the  hull  force  expression  (Eq.12)  yields 


I]  <vV 


dS 

+ 


(14) 


since  t  and  S  are  independent. 
S  i  nee 


cf,Q+<Pp  =  0  on  Y=0 

the  region  of  integration  can  be  extended  over  the  waterplane  area  Sj  . 
Thus 


-pf  ft  -fl 

S+Si 


dX. 

i 

on 


dS 

n=0+ 


The  above  form  is  suggestive  of  Green's  second  identity  since  both  S^+Cp 
and  X.  are  harmonic  functions  in  the  inner  region  and 
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and  also 


Thus 


9p+  (on  outside)  =  9p  (on  inside) 

X  =  V 


and  n  =  -n 
+ 


(1)  5  pp  3  .  . 

Fi  =  +Pf  37  II  3^ 

S+Si 


X  J  ds 
n=0_ 


(15) 


The  boundary  condition  given  by  Eq.(llb)  applied  on  the  exterior  of  the 
surface  S  yields 

3 


37  (VV|+  ~  ° 


(16) 


Furthermore,  the  properties  of  distributions  of  sources  on  curved  boundaries 

*  ,J0 

yield 


3n 

39 


D 


(17a) 


09  1 


(17b) 


where  S+Si -Q  indicates  the  region  excluding  the  location  of  the  source 
at  X;  and  a  =  0  on  ^  (waterplane) . 

From  these  relations  (17a, 17b) 


=  4tt  a 


S90 

S:f,o  ! 

37T 

+  "  <5rT  1 

or 

% 

_  d*D 

3FT 

_  ~  3n 

S  i  nee 

,!!e| 

SIT 

_  on  1 

Eq . ( 1 8)  can  be  written  as 

a(vV  1 

-  4tio 


(18) 


TFT 


rD r^P' 
Tn 


-  4rro 
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which  upon  substitution  into  Eq .(15)  yields 

Fi  '  J|axids  *  “f  ST  JJ  k  (VV 

Here  use  was  made  of  the  fact  that 

Xj=  0  on  waterplane  Y  =  0 

and 


dS 


3n  (cpD+CpP) 


=  0  on  S  (wetted  surface) 


(19) 


In  addition,  on  the  waterplane  surface  Y  =  0 


8_ 

Bn  _ 


B 

By 


and  hence 


(VV 


57  (VV 


Thus  Eq.(19)  giving  the  force  component  F;  becomes 

F,  ■  -4"pf  It  JIS  °XidS  -  "f  It  JI  xt  ly  <VVlv,0ds  (20) 

It  should  be  noted  that  the  second  term  on  the  right-hand  side  of 

Eq.(20)  is  zero  for  the  vertical  force  since  Xg=Y=0  on  ^  .  Consequently, 

when  there  is  a  free  surface  this  term  is  present  only  if  the  transverse 

and  longitudinal  forces  are  required.  Calculations  have  shown  that  the 

primary  contribution  to  the  total  unsteady  surface  force  is  the  temporal 

part  given  by  Eq.(20).  The  convective  term  described  next  may  be  neglected 

4 

as  was  found  by  Vorus,  and  by  the  present  numerical  calculations. 

The  convective  term  of  the  extended  Lagally  Theorem  is  given  by 


(2)  "  &  — 1  . 

Fi  =  'PJI  LVsj  0x7  ^VD+V  J  STT  dS 

^  J 


(21) 


which  is  the  second  term  in  Eq . ( 1 2) .  In  s teady-s tate,  propel ler/hu 1 1 
interaction  problems,  the  convective  term  described  by  Eq .(21)  is  the  only 
contribution  to  the  surface  force  on  the  hull.  Modeling  the  hull  by  a 
similar  type  of  surface  distribution  of  singularities  given  by  Eq .(13) > 
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Cox,  et.al., 


demonstrated  that  the  convective  term  can  be  expressed  as 

a<vV 

-  -4n?f  Jj  uPi  °ods  +  Pf  n  oxf  oy - ds*  <22) 


where  =  the  steady  state  source  strength  distribution 

and  Up.  =  propeller-induced  velocity  in  the  i-direction  at  any  point 
1  on  the  hull  surface 

(See  Appendix  A  of  Reference  11.) 


Consequently,  the  total  force  on  the  body  is  given  by  the  extended 
Lagally  Theorem  as 


F. 


+ 


(23) 


If  harmonic  time  dependence  is  assumed,  then  for  the  unsteady  hull 
force,  Eq.(23)  is  written  as 


*(*N) 


F.  =  -4rrPf(i  4Nuu)  5  X.ds  -  JJ  Up.  o^dS 


acp  d($  -kp p) 

-pf  II  L'^xi  +w:\  — 

Si  ' 


dS  (24) 
Y-0 


where  the  circumflex  *  marks  the  part  of  any  quantity  q  which  is  time 
dependent,  i.e.,  setting 

q(X,Y,Z  ;t)  =  q(X,Y,Z)e'  £NJJt 

The  frequency  selected  is  the  frequency  at  which  the  propeller  velocity 
field  is  radiated.  It  is  Eq.(24)  which  is  programmed  to  compute  the  hull 
surface  forces  according  to  the  procedure  developed  in  this  study. 
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POTENTIAL  FLOW  ABOUT  THE  SHIP  HULL 

The  problem  considered  in  this  section  is  that  of  the  flow  of  an 

ideal,  inviscid  fluid  about  an  arbitrarily  shaped,  three-dimensional  ship 

hull  modeled  by  a  surface  distribution  of  sources.  The  method  developed  and 

2 

programmed  by  Hess  and  Smith  is  employed  to  solve  the  Fredholm  integral 
equation  of  the  second  kind  generated  by  the  condition  that  the  hull  is  a 
stream  surface.  The  source  distribution  is  discretized  into  plane  quadri¬ 
lateral  elements  of  constant  density  and  the  simultaneous  algebraic  equations 
resulting  from  the  discretization  of  the  integral  equations  are  solved  by  a 
Gauss-Seidel  iteration  scheme.  The  numerical  analysis  procedure  and  theoret¬ 
ical  model  employed  by  Hess  and  Smith  are  adequately  described  in  Reference  2. 
A  synopsis  of  the  theory  and  the  resulting  boundary-value  problem  is  given 
below : 

The  ship  hull  surface  denoted  by  S  is  given,  in  Cartesian  coordinates, 
by  the  equation 

F(X,Y,Z)  -  0  (25) 

where  X,Y,Z  are  defined  in  Figure  1.  The  hull  is  subjected  to  the  free 
stream  and  the  propeller-induced  velocity  fields.  The  steady-state  propel¬ 
ler  inductions  are  ignored  in  the  present  investigation;  only  the  unsteady 
(blade  frequency)  hull  surface  forces  are  sought  of  which  the  convective 
contribution  can  be  shown  to  be  small.  The  prope I  1 er- i nduced ,  blade  fre¬ 
quency  velocity  field  is  expressed  in  complex  form  (as  described  in  the  next 
section)  and  has  cosine  (real)  and  sine  (imaginary)  components.  The  real 
and  imaginary  inductions  are  considered  separately  as  two  input  flows  in  the 
Hess-Smi th  method.  Therefore,  the  three  onset  flows  about  the  ship  hull 
which  are  considered  as  input  to  the  Hess-Smith  computer  code  are:  (a)  A 

a 

uniform  onset  flow  V  =  Ui ,  (b)  Rea]  part  of  the  propeller-induced  velocity 

R~oo 

I 

neia  Vp  ,  (c)  Imaginary  part  of  the  propeller-induced  velocity  field,  Vp. 

The  Hess-Smith  potential  satisfies 

v'j*hs  =  0  (26) 

on  the  region  exterior  to  S  which  is  the  hull  boundary  plus  its  image.  In 
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addi tion, 

|i  -  0  on  S  (27) 

where  n  is  the  outward  directed  normal  to  the  hull  pointing  into  the  fluid 
outside  the  ship  hull.  If  a  free  surface  exists  (denoted  Si),  then  the 
X,Z, -plane  is  taken  to  be  the  waterplane.  It  is  then  considered  to  be  a 
plane  of  geometric  symmetry  with  a  Dirichlet  or  a  Neumann  boundary  condition 
specified  on  the  free  surface  if  either  a  high  Froude  number  or  a  low  Froude 
number  approximation,  respectively,  is  desired. 

The  total  potential  for  steady  motion  may  be  expressed  by 

$HS  =  UX  +  {PS(X,Y,Z)  (28) 

where  9<.  is  the  response  of  the  hull  to  UX  flow.  It  satisfies 

V3<ps  =0  (29) 

exterior  to  S,  and 


on  S 


(30) 


where  n  =  ±[vF/| VF 1 ]p_g  .  The  sign  of  n  is  such  that  it  points  into  the 
fluid  exterior  to  S.  Also,  it  is  required  that 

9S  -  0  as  X2  +  Y2  +  Z2  -  *  (31) 

For  a  surface  distribution  of  sources,  the  diffraction  potential  is 
given  by 

9S(X,Y,Z)  =  -  JJ  ds  (32) 

x 

where  R  °  i(X-§)2  +  (V-T))2  +  (Z-£)2|2  and  S  includes  the  image  surface. 
The  form  of  9  in  Eq.(32)  automatically  satisfies  Eqs.(29)  and  (31)  for 
any  function  a  .  In  addition,  the  function  a  must  be  determined  so  that 
the  potential  9^  satisfies  the  boundary  condition  on  S  given  by  Eq.(30). 
It  should  be  noted  that  the  surface  S  is  made  up  of  the  wetted  surface  and 
that  of  the  image  system  (positive  geometric  image)  of  the  hull,  so  that  a 
closed  body  in  an  infinite  fluid  is  generated,  a  situation  which  has  been 
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shown'*  to  be  sufficient  for  the  evaluation  of  the  vertical  force.  Applying 
Eq.(30)  requires  the  evaluation  of  the  normal  derivative  of  Eq.(32)  at  a 
point  (X,Y,Z)  on  S.  As  the  surface  S  is  approached,  the  derivative  of 
the  integral  becomes  singular  and  its  principal  part  must  be  extracted 
(see  Kellogg'^).  This  manipulation  leads  to  the  following  integral  equation 
for  the  surface  source  density  distribution: 

-n  •  Ui  =  2ttct(p)  -  JJ  a(Q)  ^  dS  (33) 

where  a  is  the  source  density,  n  is  the  normal  vector  at  P,  and  R  is 
the  Descartes  distance  between  P  (control  point)  and  Q.  (loading  point). 

This  equation  is  a  Fredholm  integral  equation  of  the  second  kind  for  the 
unknown  source  densities  a,  which  are  determined  for  the  known  onset  flow 
on  the  hull  surface  specified  on  the  left-hand  side.  The  velocity  Ui  is 
the  velocity  at  P  due  to  the  undisturbed  stream  at  upstream  infinity.  Any 
flow  field  imposed  on  the  locus  of  the  hull  (in  the  absence  of  the  hull) 
may  be  considered.  Here  two  other  fields  are  considered,  namely,  the  cwo 
components,  real  and  imaginary,  of  the  velocity  field  Induced  by  the  pro¬ 
peller.  For  each  given  onset  flow,  the  source  densities  are  computed  by 
sol vi ng  Eq. (33)  . 

In  the  present  investigation  the  coded  numerical  procedure  developed 
2 

by  Hess  and  Smith  is  applied  to  solve  Eq . (33) .  The  three  onset  flows  con¬ 
sidered  produce  three  source  densities,  namely,  oq  due  to  the  uniform  steady 
motion  of  the  hull  and  its  positive  image  moving  in  an  infinite  fluid,  a 

K 

due  to  the  real  part  of  the  prope 1 1 er- i nduced  onset  flow,  and  due  to  the 
imaginary  part  of  the  prope 1 1 er- i nduced  onset  flow.  The  blade  frequency 
radiated  velocities  induced  by  the  propeller  are  the  velocities  imposed 
here . 
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VELOCITY  FIELD  INDUCED  BY  THE  PROPELLER 

Considerable  effort  has  been  devoted  to  the  development  of  unsteady 
lifting  surface  theory  for  propellers  operating  in  a  spatially  nonuniform 
hull  wake.  The  theoretical  analysis  and  numerical  techniques  developed  at 
Davidson  Laboratory  are  based  on  modeling  the  propeller  blades  with  sources 
and  pressure  dipole  singularities.  The  computer  code  permits  the  calcula¬ 
tion  of  unsteady  blade  pressure  distributions,  forces,  and  moments  for  a 
given  propeller  geometry  —  number  of  blades,  pitch,  chordlength,  camber, 
thickness,  skew,  and  rake  —  and  a  given  inflow  field.  Details  of  the 
theory,  analysis  and  computer  program  have  been  reported  elsewhere  (see 

c  ^ 

Tsakonas,  et.al.  ’  ).  The  theory  has  been  extended  to  compute  the  field 

g  Q  12 

point  velocities.  ’  Both  the  loading  and  velocity  field  programs  have 
been  incorporated  into  the  prope 1 1 er/hu 1  1  interaction  computer  code  des¬ 
cribed  here.  The  steady  and  unsteady  blade  loadings  and  their  concomitant 
field  point  inductions  are  the  responses  of  the  propeller  blades  to  the  hull 
wake. 

The  hull  wake  longitudinal  and  tangential  components  as  measured  in 
the  propeller  plane  are  resolved  by  Fourier  analysis  into  a  circumferential 
mean  component  and  higher  harmonics  in  shaft  rotation  frequency.  The  higher 
harmonics  which  are  oscillatory  in  time  with  respect  to  the  propeller  give 
rise  to  unsteady  loadings  on  the  blade  in  a  manner  analogous  to  an  airfoil 
encountering  a  sinusoidal  gust.  Within  the  linear  approximation,  each  pro¬ 
duces  a  component  of  unsteady  blade  loading  at  the  corresponding  frequency. 

Upon  summing  over  all  the  blades,  only  certain  harmonics  contribute  to  the 
net  exciting  forces  and  moments  acting  on  the  shaft,  namely,  £N  for  the 
thrust  and  torque  and  £N±1  for  the  side  bearing,  forces  and  bending  moment, 
where  l  is  any  integer  and  N  is  the  number  of  blades.  However,  _a_l_l_  the 
harmonics  of  loading  contribute  to  the  forces  on  an  individual  blade  and 
to  the  blade  rate(and  integer  multiples  thereof)  of  radiated  velocity  and 
radiated  pressure  fields  of  the  propeller  (e.g.,  at  points  on  the  hull 
surface)  . 

A  brief  description  of  the  propeller  theory  applied  to  compute  the  field 
point  velocities  (which  includes  the  velocities  induced  by  the  blade  thickness 
modeled  by  source  s i ngu 1 ar i t i es)  is  presented  next; for  further  details  the  reader 
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may  review  References  5,  6,  8,  and  9.  The  loading  program  developed  at  the 
Davidson  Laboratory  is  called  PPEXACT.  Prior  to  calculating  the  hull  vibra¬ 
tory  surface  forces  due  to  the  propeller  by  exercising  the  propel ler/hul 1 
interaction  program  PIHF,  the  program  PPEXACT  (or  the  like)  must  be  exercised 
to  determine  the  blade  loadings  for  all  frequencies  up  to  one  plus  the  blade 
frequency,  and  also  to  determine  the  unsteady  shaft  forces  and  moments.  The 
blade  frequency  hull  forces  are  computed  by  PIHF  which  uses  the  blade  load¬ 
ings  as  input  to  compute  first  the  veloc i ty  field  produced  by  the  propeller. 

The  algorithm  in  the  PIHF  code,  which  utilizes  the  field  point  velocity 
program,  the  Hess-Smith  program  and  the  extended  Lagally  Theorem  code,  is 
described  in  the  next  section. 


The  propeller  theoretical  analysis  uses  the  acceleration  potential 
method  to  formulate  a  linearized,  unsteady  lifting-surface  theory  for  a  marine 
propeller  operating  in  a  spatially  varying,  unidirectional  mean  flow  which  is 
aligned  with  the  shaft  axis.  Potential  flow  of  an  ideal,  inviscid,  incompres¬ 
sible,  homogeneous  fluid  is  assumed.  The  steady  forces  and  moments  produced 
by  a  propeller  operating  in  a  real  fluid  are  due  primarily  to  potential  flow 
effects;  however,  since  the  viscous  drag  is  significant  in  predicting  steady 
performance,  it  is  accounted  for  approximately  in  the  PPEXACT  program.  The 
unsteady  forces  and  moments  are  considered  to  be  inviscid  in  nature.  The 
vorticity  in  the  oncoming  flow  is  ignored;  however,  the  variations  in  the  in¬ 
flow  velocity  are  considered  parametrically  as  is  usually  done  in  theories  for 
marine  propellers.  The  blades  are  assumed  to  lie  on  helicoidal  surfaces. 
Therefore,  the  linearized  equation  of  motion  for  unsteady  flow,  referred  to 
a  nonrotating  (fixed)  cylindrical  coordinate  system  ( x , r , c)  centered  at  the 
propeller  axis  (see  Figure  2)  may  be  written  as 


di 

at 


+ 


(34) 


where  |  is  the  so-called  acceleration  potential  and  p  is  the  local  pertur¬ 
bation  pressure.  Hence, 


V'p  =  0  (35a) 

p—0  as  x  —  -  ®  (35b) 


The  solution  to  Eq.(34)  subject  to  (35b)  can  be  shown  to  be 


f(x,y,z;t)  =  -  —  j'p(T,y,z;t  +^)d" 

I  -00 
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where  p  is  the  disturbance  pressure  due  to  the  propeller  loading,  and 
$(-»,y,z)=0  was  used. 

The  propeller  blades  are  modeled  by  pressure  dipoles  (as  far  as  their 
loading  is  concerned),  and  the  corresponding  pressure  field  is  given  by 

p =  -  h  n  ap  h  (i)ds  (37) 

s 

After  an  harmonic  analysis  of  Ap,  Equation  (37)  may  be  written  as 

f-k  JI£»eW<5,P,0)e'Wt§^(i)dS  (38) 

where  R  is  the  Descartes  distance  between  the  control  point  (x,r,cp)  (i.e., 
fixed  point  in  space)  and  the  loading  point  (§,p,0)  rotating  with  angular 
velocity  -0  and  lying  on  a  helicoidal  surface. 

Equations  (36)  and  (38)  are  combined  and  summed  over  the  number  of 
blades  to  yield  the  equation  for  the  velocity  potential  at  a  point  (x,r,cp) 
due  to  the  blade  loading  of  an  N-bladed  propeller,  namely, 


$(x,r,cp;q)=  + 


1 

4TTpfU 


N  » 


AP(X)(S,P,0o)e 


x(cit-en) 


X'  i?,a(r-x) 
,  e 


(39) 


where  q  and  \  are  positive  integers  defining  the  desired  frequency  of 
the  velocity  field  and  the  appropriate  frequency  of  the  propeller  loading, 
respectively.  Note  that  q=0  corresponds  to  the  steady-state  velocity  field, 
and  q=N  corresponds  to  the  blade  frequency  induced  velocity  field.  The 
latter  is  of  interest  here.  The  other  parameters  are  defined  as  follows: 

x,r,cp  =  cylindrical  coordinates  of  point  in  space  referred  to 
axes  with  origin  at  the  propeller  hub 

§,p,9o=  cylindrical  coordinates  of  loading  point  on  propeller 
U  =  free  stream  velocity,  fps 
Pj.=  fluid  density,  slugs/cu  ft 
-Li  =  angular  velocity  of  the  propeller,  rad/sec 
a  =  w/u,  1/ft 

AP^(f,p,0  )=  propeller  loading,  or  pressure  jump  across  the  propeller 
blade,  psf 
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R  =  |(T-§)2  +  r  +  p3  -  2 rp  cos[eo-Qt  +  en  -  a(T-x)-cp]  i2,  ft 

- —  =  normal  derivative  on  the  helicoidal  surface  at  the  loading 
°n  point  (|,p,0o)  on  the  propeller 

-  2tt 

0n  =  —  (n-1)  n  =  1,2,  ...  N 
S  =  blade  surface,  sq  ft 

The  normal  derivative  to  the  helicoidal  surface  (specified  by  §  =  0o/a)  is 


d 

377 


i  5  1  d 

la  If  --J6Q-J 


(1  +  a2  p2 ) 2 
Use  is  made  of  the  following: 

a)  Expansion  scheme  for  the  reciprocal  of  the  Descartes  distance  R 


1 

R 


i  Z7  eIm(e0-at+®n-a(T-x)"cP) 

n  m=  -  °° 


Jlm(lkl  P)  Km(  l  k  I  r)  e 1  ^  ^  kdk 


(when  p  <  r,  otherewise  p  and  r  are  interchanged  in  the  modified  Bessel 
funct i ons) . 

b)  Summation  over  the  propeller  blades 

N  -i(\-m)0n  ,  N  when  7-m  =  £N,  Z  =  0,-1, -2,  ... 

Li  e  =  j 

n=l  t  0  when  \-m  /  XN 

and  from  the  time-dependent  factor 
X  -  m  =  q  =  £N 


c)  T  ransfo rmations 

L(x)(p,ea)  =  APU)(§,p,e0)-peb,  ib/ft 

where  0^  is  the  projected  semichord  at  each  spanwise  location  p,  in  radians. 
0o  =  °  "  Gb  cos  °a  =  a? 

where  a  is  angular  position  of  the  midchord  line  of  the  projected  blade 
from  the  generator  line  through  the  hub,  or  skewness,  and  bQ,  is  angular 
chordwise  location  of  the  loading  point. 
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d)  Approximation  of  the  loading  function  (P>®q.)  ‘n  the  chordwise 
direction  by  n  modes  whose  selection  is  dictated  by  the  pressure 
distribution  on  a  foil  in  two-dimensional  flow 


The  velocity  potential  can  then  be  written  as 

$ 


CO  CO 


<x,r, <?;,)-  JE  £  L<X’">(P)eit'Q,C'"'n>dp 
p  n= I  A=o 


m 


where  the  function  K^m,n^,  which  is  known  as  the  propagation  or  influence 
function,  is  given  by 

- (m, n)  iNe'm^CC^  /  iiN(a-ax),'m  2,„i 

K  "  '  -frrpfUa-  -  le  A7  "  a  iH) 

•  1  (a|iN|p)Km(a|XN|r)A(n)(\0b)  -  ~  j  e‘ k(x_CT/a) 


(m/p2 +ak) 
k+a£N 


lm(tklP)Km(|k|r)-AlnJ^m-  -;6bjdk/  (4l) 


for  p  <  r  (for  p  >  r,  p  and  r  are  interchanged  in  the  modified  Bessel  func¬ 
tions).  The  symbol  (x)  denotes  the  results  of  the  chordwise  (9^)  in¬ 

tegration  which  differs  with  the  assumed  chordwise  modes. 

The  following  chordwise  mode  shapes  have  been  used  to  model  blade 
pressure  distributions  when  A=0  ; 

(a)  The  so-called  "roof-top"  mode  appropriate  to  NACA-a-mean 1 i ne 
cross-sections  when  the  propeller  operates  at  the  design  advance  ratio.  The 
chordwise  loading  is  constant  from  x/c  =0  at  the  leading  edge  to  x/c  =  a 
and  then  decreases  linearly  to  zero  at  x/c  =  1,  the  trailing  edge. 

(b)  The  first  term  of  the  Birnbaum  series  (the  cotan  0^/2  term)  for 
the  additional  loading  due  to  difference  between  off-design  and  design 
advance  ratio  when  the  propeller  operates  at  off-design  conditions. 

(c)  The  complete  Birnbaum  series  when  the  design  operating  condi¬ 
tions  are  not  known. 

When  X/0  ,  the  blade  loading  distribution  is  represented  by  the  complete 
B i rnbaum  se r i es . 
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The  velocities  in  the  x,  y,  and  z  directions  of  the  Cartesian  co¬ 
ordinate  system  with  origin  at  the  propeller  center  at  a  point  x,r,cp  due  to 
the  propeller  operating  in  spatially  varying  inflow  will  be  given  by  the 
appropriate  components  of  the  gradient  of  the  velocity  potential  $. 

The  axial  component  in  nondimens ional  form  is 


U  U 


and  the  y  and  z  components  are  derived  from 

V  §  §  ,4 

_X  =  =  _  _L  sincp  .  1  J£  coscp 

U  U  u  y  r  u  y 


V  $  § 

z  z  r 

“  =  IT  =  IT  coscp 


r  U 


s  i  ntp 


Reference  8  shows  that  for  the  unsteady  blade-frequency  case, 
q=N,  \l\  =  1,  which  is  of  interest  ?n  the  present  study, 

u  u 


v[N)  *  (N) 


n  max 

”  I  -s, 

p  n=l 

00 

£ 

\®0 

jLU,n) 

l 

+  conj . 

[L(x,S)((j)^«.S)]}da 

n  max 

00 

-  I  s 

E 

f  (X, in) 

r 

P  n=l 

\=0 

+  conj . 

[La.5>(il)R(WN,s)’}  dp 

viN)  V10 


n  max  « 


7j  S2  (p) -(X-N,n) 

p  n=l  \=0  ^ 


+  conj.[LU’n)(p)K^+N,n)j}  dp  (42) 


The  propeller-induced  velocity  field  is  thus  made  up  of  a  large  number  of 
combinations  of  the  frequency  constituents  of  the  loading  function  with 
those  of  the  propagation  function. 
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The  loading  functions  L^’n^(p)  are  determined  through  the  PPEXACT 
program  which  solves  the  integral  equation  relating  the  unknown  blade  load¬ 
ing  distribution  with  the  known  onset  flow  normal  to  the  blade  surfaces. 

The  derivatives  of  the  propagation  functions  K^m’^  present  a  complex 
problem.  Beside  their  oscillatory  character  in  terms  of  the  spatial  or¬ 
dinates  which  presents  the  problem  of  selecting  the  proper  spacing  or  mesh 
size,  there  is  the  problem  of  truncation  of  the  range  of  integration. 


In  deriving  the  velocity  potential  at  a  field  point  x,r,cp  due  to 

blade  thickness  of  an  N-bladed  propeller,  the  thickness  distribution  of 

a  blade  section  is  represented  by  a  source-sink  distribution  smeared 

9 

over  the  blade  section  chord.  This  velocity  potential  is  given 
by 


*T(x,r,cp;t)  = 


1 

2+ir 


N 

E 

n=l 


b  M(|, p, 0Q)  VJ 


i  J 


-0 


P  P 


+a2  p3 
ap 


pdpdB 


(43) 


The  source  strength  H ( § , p , ©0)  is  determined  in  accordance  with  the  "thin 
body"  approach  as 

df(5,P,ej 

M(S,P,e0)  =  2u  -  -^r  - 


where  f(|,p,0o)  is  the  thickness  distribution  over  the  blade  section  at 
the  radial  distance  p  .  The  Descartes  distance  from  a  point  on  the  propel¬ 
ler  to  the  field  point  is 

X 

Rt  =  i(x-§)S  +  rs  +  p2  -  2rp  cos(eo-Qt-cp+en) \ 2 

It  is  known  that  the  detailed  geometry  of  the  blade  has  little  or  no 
effect  on  the  velocity  or  pressure  field  around  an  operating  propeller.  This 
is  in  contrast  to  the  case  of  the  pressure  distribution  on  the  propeller 
blades,  where  accurate  and  detailed  information  about  the  propeller  blade 
section  geometry  including  the  thickness  distribution  is  important. 

Assuming  that  the  blade  section  is  bounded  by  two  circular  arcs  which 
is  equivalent  to  saying  that  the  section  is  of  the  lenticular  form,  then 
the  thickness  distribution  can  well  be  represented  by 


li 
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f(P.0a)  =  slr>2Qoi 


t  (p)  0J”  , - 

si"% 


where  t(p)  is  the  maximum  thickness  and  t  /c  the  maximum  thickness  ratio 
v/  o 

to  expanded  chord. 

Utilizing  the  expression  for  the  thickness  distribution  together  with  the 

known  expansion  of  the  reciprocal  of  the  Descartes  distance,  the  velocity 

g 

potential  for  the  propeller  blade  thickness  becomes 


$  (x,r,cp;t)  £  -iAN(cp+Ot)  f  eiiNgP 

T  n  2  i—  ji  0? 


•  I  'ifjCk'PjKxi N(lk(r)eik(x'aP/a) 


{si  n[  (k-aJlN)  Q^/a]-  (k-a^N)  ( 0  P/a)  cost  (k-a£N)  0^/a]  (dkdp 

(k-ai.N)3 


m 


The  velocity  components  along  the  axes  of  the  Cartesian  coordinate  system 
due  to  thickness  can  be  determined  by  taking  the  appropriate  derivatives 
and  these  velocities  are  then  combined  with  those  due  to  the  loading  to 
obtain  the  total  propeller-induced  velocity  field. 

A  computer  program  adapted  to  a  CDC  6600  or  7600, or  Cyber  176  high¬ 
speed  digital  computer  has  been  developed  for  the  evaluation  of  the 
propeller-induced  velocity  field. 

The  analyses  of  References  8  and  9  and  the  corresponding  computer 
program  have  been  modified  to  furnish  the  contribution  of  the  propeller  image 
system  to  this  velocity  field  at  high  Froude  number.  The  use  of  the  negative 
of  the  propeller  potential  located  at  the  reflected  position  above  the  free 
surface  satisfies  the  high  Froude  number  condition  on  that  surface.  For 
points  on  the  free  surface  this  condition  yields  zero  perturbation  pressure 
or  equivalently  doubles  the  vertical  component  of  the  velocity  induced  by 
the  direct  propeller.  The  computer  program  automatically  gives  the  propeller- 
induced  velocity  field  generated  from  the  direct  and  image  systems  together  or 
separately,  when  the  high  Froude  number  condition  on  the  free  surface  is  to 
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be  satisfied. 

The  program  requires  the  coordinates  of  the  field  points  to  be  given 
in  cylindrical  coordinates  x,r>9: 

a)  For  the  direct  system 


where  the  angular  ordinate  9^  is  measured  in  the  counter-clockwise  direc¬ 
tion  from  the  vertical  (12  o'clock)  position  looking  forward,  and  x,y,z 
are  the  rectangular  coordinates  of  the  axes  fixed  in  the  propeller. 


b)  For  the  image  system 


r(  =  /y2  +(2d-z)2 

^1  =  tan_1  (ife) 


and  is  measured  clockwise  from  the  6  o'clock  position  looking  forward; 

d  =  distance  of  the  propeller  axis  from  the  free  surface. 

The  resolutions  in  the  y  and  z  directions  are  given  as  follows: 

a)  For  the  di rect  system 


b) 


^VD  _  Vro  .  V 
U  “  U  Sl%  +  u 


%  Vr 


U  U  ‘''~TD  U 
For  the  image  system 

,  _j_  .  V' 

—  -I-  c  1  r»rn  J> 

u 


VV|  V_I±  . 

- u  5,n*i  +  U 


V 


Z|  VZ|  V9; 

u  =  u  cos^i+  u 


cos9D 

sin9D 

cos9j 

sin9, 
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In  the  present  problem,  the  blade  frequency  velocities  due  to  the 
propeller  and  its  image  are  computed  at  selected  hull  points  where  their 
contributions  are  considered  to  be  significant.  The  velocities  imposed 
by  the  propellers  on  the  hull  in  the  vicinity  of  the  propellers  are  com¬ 
puted  first.  Then  the  range  of  computed  points  is  extended  until  it  is 
observed  that  the  additions  to  the  force  are  negligible  from  the  most 
distant  panels  considered.  All  the  quadrilateral  elements  representing 
the  hull  that  are  outside  the  direct  range  of  propeller  influence  can 
therefore  be  set  to  zero. 

The  real  and  imaginary  parts  of  the  prope 1 1 er- i nduced  velocities 
calculated  by  the  procedure  described  above  yield  the  two  sets  of  boundary 
conditions  ("real"  and  "imaginary"  onset  flows)  applied  to  the  ship  by 
input  to  the  Hess-Smi th  program,  which  is  part  of  PIHF. 

4 

Vorus  has  noted  that  the  calculation  of  just  the  vertical  force 
can  be  accomplished  by  treating  the  hull  and  its  image  as  a  completely 
submerged  body  in  the  presence  of  the  propeller  alone,  this  be:ng  equiva¬ 
lent  to  the  force  generated  on  the  wetted  portion  of  the  hull  in  the 
presence  of  both  its  geometric  image  (containing  sources  of  opposite  sign) 
and  the  propeller  plus  its  negative  image.  As  this  observation  obviates 
the  need  for  the  inductions  from  the  image  propeller,  we  have  calculated 
the  vertical  force  as  though  the  body  were  completely  submerged  in  the 
presence  of  only  the  propeller.  The  first  term  of  Eq .(12)  applies  with  S 
representing  the  double  body  surface. 
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CORRELATION  OF  THEORY  WITH  EXPERIMENTS 


By  combining  programs  for  computing  the  unsteady  blade  loadings  and 
the  field  point  inductions  with  a  program  for  the  evaluation  of  hydrodynamic 
reactions  of  a  hull  to  those  inductions,  a  theoretical  procedure  has  been 
developed  and  a  new  program  has  been  established  and  adapted  to  a  CDC-6600 
series  or  Cyber  176  high-speed  digital  computer  which  furnishes  (a)  the 
propel  1 er- i nduced  velocities  at  points  on  a  given  hull,  (b)  the  hydrody¬ 
namic  reactions  of  the  hull  to  those  inductions,  and  (c)  the  unsteady  forces 
generated  on  the  hull  through  the  bearings  and  the  hull  plating  (by  the  fan¬ 
ning  action  of  the  blades).  The  accuracy  and  usefulness  of  the  computational 
procedure  is  investigated  below  by  comparing  theoretical  predictions  with  ex¬ 
perimental  measurements  obtained  in  two  widely  different  cases. 

The  first  set  of  model  data  was  provided  by  measurements  made  at  CWT 
Naval  Ship  Research  &  Development  Center  on  a  fully  submerged  spheroidal  head 
form  cantilevered  from  a  force  balance  as  illustrated  schematically  in  Fig¬ 
ure  2.  The  force  measurements  were  made  with  a  three-bladed  propeller 
NSRDC-4118  located  16  inches  aft  of  the  nose  of  the  body  with  a  nominal  tip 
clearance  of  3  inches  operating  at  advance  ratio  J  =  0.833  (design  speed). 

The  second  case  selected  was  a  Series  60  hull  model  (v-stern)  driven 

by  the  four-bladed  propeller  NSRDC  3379  which  was  tested  by  the  late  Profes- 

13 

sor  Emeritus  F.M.  Lewis  at  MIT.  A  schematic  arrangement  of  the  propeller- 
hull  configuration  and  the  coordinate  axis  system  are  presented  in  Figure  1. 

To  carry  out  the  numerical  solution  of  the  Fredholm  integral  equation 
of  the  second  kind  given  by  Eq  .(33) >  the  hull  surfaces  are  subdivided  into 
quadrilateral  elements  as  shown  in  Figures  3  and  4.  The  quadrilateral  size 
is  crucial  in  calculating  the  hydrodynamic  forces  and  moments,  but  it  is 
rather  difficult  to  select  the  proper  dimensions.  These  are  dictated  by  the 
propeller-induced  velocity  which  is  of  oscillatory  character  as  shown  in  Fig. 5. 
The  patch  dimensions  should  be  selected  so  as  to  fit  four  quadrilaterals 
into  one  cycle  of  spatial  oscillation  wherever  such  variations  are  present. 

It  has  been  suggested  by  Dr.  Breslin  that  rather  than  calculate  the 
prope 1 1 e r- i nduced  velocity  components  in  advance,  an  approximate  idea  of  the 
periods  of  oscillation  could  be  obtained  from  the  exponential  factors 
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1 

exp(ilNax)  and  exp(  i  ( X.-£N) cp)  .  Then  for  fixed  radial  clearance  and  angular 
position  cp,  the  longitudinal  spacing  in  the  neighborhood  of  the  propeller 
would  be  governed  by  Nax=2n/4  or  x=2n/4Na  .  For  fixed  longitudinal 
position  the  transverse  spacing  would  be  governed  by  | \-iN|9=2n/4  which 
for  .£=1  and  \=0  (uniform  flow)  yields  cp  =  2ti/4n. 

t 

In  the  case  of  the  body  of  revolution  the  mesh  size  should  be  governed 
by  x~0.l4  propeller  radius  and  cp~  0.  5  so  that  the  quadrilateral  sizes  se-  | 

lected  (see  Fig. 3)  are  adequate.  In  the  case  of  the  surface  ship  the 
longitudinal  spacing  should  be  x~0.08  propeller  radius  and  the  transverse 
spacing  governed  by  cp~0.4.  The  longitudinal  spacing  shown  in  Figure  4  is 
inadequate  but  the  transverse  spacing  may  be  considered  sufficient.  To  im¬ 
pose  the  required  longitudinal  spacing  would  increase  the  already  large 
number  of  quadrilaterals  considerably  and,  hence,  the  number  of  field  points 
at  which  the  three  velocity  components  are  to  be  evaluated  with  all  the  com¬ 
binations  of  the  frequencies  of  the  loading  and  the  propagation  functions. 

This  is  impractical. 

A  study  has  been  initiated  at  Davidson  Laboratory  to  replace  the  pres¬ 
ent  method,  which  represents  the  hull  by  a  source  distribution  and  requires 
computation  of  three  velocity  components  at  each  field  point, by  a  procedure 
which  represents  the  hull  by  a  doublet  distribution  and  requires  the  compu¬ 
tation  of  only  one  function,  the  propeller  velocity  potential.  This  method 
permits  decreasing  the  mesh  size  of  the  quadrilaterals  and  increasing  the 
number  while  still  reducing  computer  time,  and  is  showing  considerable  improve¬ 
ment.  However,  for  this  report  the  calculations  have  been  performed  using  the 
quadrilateral  arrangements  shown  in  Figures  3  and  4. 

CASE  1 — The  surface  of  the  spheroidal  body  is  divided  into  154  quadri¬ 
laterals  with  mesh  size  increasing  with  distance  from  the  propeller.  Each  of 
these  elements  (see  Figure  3)  is  a  uniform  source  patch.  Panels  93  through 

100  were  used  to  close  the  body.  The  geometry  of  these  elements,  together  I 

with  the  normal  velocity  (due  to  loading  and  blade  thickness)  induced  by  the 
propeller  operating  in  uniform  inflow,  form  the  input  to  the  generalized 
Hess-Smith  program  which  inverts  the  integral  equation  (similar  to  Eq.33) 
and  evaluates  the  source  densities  on  each  of  the  panels. 

The  particulars  of  the  propeller  are  given  in  Table  1  and  other  charac¬ 
teristics  are  given  in  Table  2. 
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TABLE  1 


PARTICULARS  OF  PROPELLER  A 118 


Propeller  Designation 

NSRDC  41 18 

Number  of  Blades 

3 

EAR 

0.6 

P/D  at  0.7 

1 .08 

Diameter,  D,  ft 

1 

RPM 

900 

Fluid  Speed,  U,  ft/sec 

12.5 

Advance  Ratio,  J  =  U/nD 

0.833 

TABLE  2 

CHARACTERISTICS  OF  PROPELLER  4118 


z/zo 

Max i mum 
camber/chord 
mx/c 

Maximum 

thickness/chord 

fo/c 

Ratio  of  leading- 
radius  to  chord 

Po 

0.25 

0.0228 

0.090 

0.00525 

0.35 

0.0231 

0.068 

0.00290 

0.45 

0.0224 

0.052 

0.00170 

0.55 

0.0212 

0.040 

0.00100 

0.65 

0.0203 

0.031 

0.00060 

0.75 

0.0198 

0.024 

0.00035 

0.85 

0.0189 

0.018 

0.00025 

0.95 

0.0174 

0.016 

0.00020 

Then  by  making  use  of  the  extended  Lagally  Theorem  (see  Eq.24),  the 
blade  frequency  force  exerted  on  the  spheroidal  head  in  the  fully-submerged 
case  will  be  given  by 

Fn  =  _4Tlpflt  II  roe’NUtdS  -  4npf  JJ  \/pe'NQto0  dS  (45) 

S  S 

where  the  first  term  represents  the  time-dependent  and  the  second  the  con¬ 
vective  term  of  the  force,  and 

r  =  position  vector  in  the  direction  of  the  desired  force 
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Vp  =  propeller-induced  velocity  vector 

ct  =  amplitude  of  the  blade  frequency  source  density 

a  =  amplitude  of  the  source  density  arising  from  the  presence  of 
the  hull  in  the  free  stream 


S  =  the  entire  wetted  surface 


In  applying  the  foregoing  formula,  it  was  found  that  the  convective 
term  is  negligibly  small  and  so  the  term  involving  the  time  derivation  can 
be  considered  the  sole  contributor  to  the  blade  frequency  hull  force.  The 
time-dependent  force  (blade  frequency  force)  in  the  i-direction  is  given 
in  the  program  by 


and 


jr.  =  -i4TTpfNaU2ro2  JJ  ?*  o*dS* 


where  n=rps 


(46) 


and  where  all  the  "starred"  quantities  are  nond i mens i ona 1 i zed  with  respect 
to  propeller  radius  rQ. 

The  theoretical  results  are  given  in  the  following  Table  3  and  the 
comparison  with  experiments  is  graphically  exhibited  in  Figure  6.  The 
graphical  comparison  shows  agreement  in  magnitude  of  theoretical  predictions 
with  corresponding  experimental  measurements  within  11%,  and  good  agreement 
in  phase  angle  9p.  9^  is  generally  defined  as  the  reference  blade  position 

when  the  force  on  the  body  is  maximum  positive  away  from  the  propeller. 


TABLE  3 

CALCULATED  BLADE  FREQUENCY  VERTICAL  FORCE  COEFFICIENT  FOR  PROPELLER  4118 
LOCATED  16-IN  AFT  OF  A  SPHEROIDAL  BODY  WITH  TIP  CLEARANCE,  c  =  3.0  INCHES 


Contributors 

F  .Vertical  Force, 
Y  lb 

Rea)  Imaginary 

CF= 

pfn  D 

9p=Phase  Angle  (Mechanical) 
in  degrees  wi th  respect  to 
Propel ler 

Coordi nates 

Th i ckness 

1.263 

-0.0133 

2.820 

-0.2° 

Load i ng 

-0.0920 

-0.2078 

0.507 

-37.96° 

Loading  + 
Thickness 

1.171 

-0.2211 

7.662 

-3.56° 
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In  this  particular  case  with  the  experiments  conducted  with  the  pro¬ 
peller  located  on  the  portside  of  the  spheroidal  body,  the  measured  and 
theoretically  defined  phase  angles  coincide. 

It  is  to  be  noted  that  the  contribution  due  to  blade  thickness 

\k 

dominates  the  force  amplitude.  This  was  also  found  by  Breslin  who  calcu¬ 
lated  the  forces  and  moments  on  infinitely  long  cylinders. 

CASE  2  —  The  Series  60,  0.6c„,  hull  model  (V-stern)  is  similarly 

-  D 

subdivided  into  quadr i la tera ] s .  There  are  202  in  all  with  the  hull  re¬ 
flection  in  the  waterplane  included.  In  some  of  the  quadrilaterals  the 
propeller-induced  velocities  are  neglected  as  being  small.  In  fact,  96  of 
the  quadrilaterals  located  at  the  forward  part  of  the  hull  and  its  reflec¬ 
tion  have  thus  been  deactivated  (i.e.,  propel  1 er- i nduced  velocities  are 
taken  to  be  zero) . 

Particulars  of  the  Propeller  NSRDC  3379  are  given  in  Table  4,  and 
other  characteristics  are  given  in  Table  5;  the  harmonic  content  of  the 
wake  behind  the  Series  60,  0.6CD  hull  with  V-shaped  sections,  is  pre- 

D 

sented  in  Table  6. 


TABLE  4 

PARTICULARS  OF  PROPELLER  3379 


Propeller  Designation 
Model  Designation 
Number  of  Blades 
EAR 

P/D  at  0.7  radius 
Diameter,  D  ,  ft 
RPM,  60n 

Model  Speed,  U  ,  ft/sec 
J  =  U/nD 


NSRDC  3379 

Series  60,  C  =0.6,  V-form 

D 

4 

0.471 
1 .02 
0.4167 
787.5 
4.56 
0.834 
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TABLE  5 


CHARACTERISTICS  OF  PROPELLER  3379 

Maximum 

Maximum 

camber/ chord 

thi ckness/chord 

a  =  skewness  in 

rad i ans 

z/z 

m  /c 

t  /c 

o 

X 

0 

0.25 

0.0254 

0.1651 

0.023 

0.35 

0.0254 

0.1303 

0.060 

0.45 

0.0263 

0.1044 

0.105 

0.55 

0.0265 

0.0846 

0.143 

0.65 

0.0257 

0.0680 

0.160 

0.75 

0.0240 

0.0543 

0.190 

0.85 

0.0207 

0.0435 

0.232 

0.95 

0.0174 

0.0334 

0.275 

TABLE  6 

HARMONIC  CONTENT 

OF  THE  WAKE"  OF  A 

SERIES  60, 

0.6  CD  M00EL  WITH  V-F0RM  SECTIONS 

H 

a)  Long! 

itudinal  Velocity 

Components 

Radius 

ao  at 

a  2  33 

a4 

a  3 

bq 

0.25 

0.522  -0.093 

-0.130  0.013 

0 

-0.005 

0 

0.35 

0.640  -0.110 

-0.190  0.027 

-0.047 

0.019 

0 

0.45 

0.736  -0.104 

-0.220  0.023 

-0.075 

0.027 

0 

0.55 

0.791  -0.142 

-0.192  -0.012 

-0.080 

0.020 

0 

0.65 

0.813  -0.152 

-0.165  -0.040 

-0.071 

0.008 

0 

0.75 

0.818  -0.158 

-0.147  -0.050 

-0.057 

-0.002 

0 

0.85 

0.828  -0.162 

-0.133  -0.057 

-0.046 

-0.012 

0 

0.95 

0.845  -0.172 

-0.122  -0.062 

-0.037 

-0.018 

0 

b)  Tangential  Velocity  Components 

Radi  us 

Aq  *i 

B  3  B3 

b4 

Bo 

0.25 

0  -0.040 

-0.040  0.037 

-0.020 

0 

0.35 

0  -0.093 

-0.035  0.030 

-0.019 

0.010 

0.45 

0  -0.128 

-0.033  0.019 

-0.017 

0.013 

0.55 

0  -0.133 

-0.034  0.005 

-0.014 

0.010 

0.65 

0  -0.130 

-0.036  -0.007 

-0.009 

0.005 

0.75 

0  -0.130 

-0.038  -0.016 

-0.008 

-0.002 

0.85 

0  -0.128 

-0.040  -0.020 

-0.007 

-0.006 

0.95 

0  -0.125 

-0.041  -0.022 

-0.007 

-0.009 

whe  re 


•'•Grant,  J.W.  and  Lin,  A.C.,  "The  Effects  of  Variations  of  Several  Parameters 
on  the  Wake  in  Way  of  the  Propeller  Plane  for  Series  60-0.60  Cg  Models"  Naval 
Ship  Res.  and  Dev.  Center  Rpt.  No. 3024,  June  1969. 


29 


R-1979 


Longitudinal  Velocity  V  /U  =  a  +  [a  cosqO  +  b  sinac)] 

x  °  q=i  q  q 

CO 

Tangential  Velocity  VT/U  =  A  +  1j  La  cosqO  +  B  sinqOj 

t  0  q=]  q  q 

Utilizing  the  extended  Lagally  Theorem  over  the  double-hull  body,  the 
vertical  component  of  the  blade  frequency  force  has  been  calculated  as  be¬ 
fore,  demonstrating  again  that  the  convective  term  is  small.  The  theoretical 
results  are  given  in  Table  7. 

TABLE  7 

CALCULATED  BLADE  FREQUENCY  VERTICAL  FORCE  COEFFICIENT 
ON  A  SERIES  60,  0.6  Cg  HULL  (V-FORM)  WITH  THE  4-BLADED  PROPELLER  3379 
AT  AXIAL  CLEARANCE,  x  =  0.25D 


Contributors 

F  =Vertical  Force, 

Gp=phase  Angle  (Mechanical) 
in  degrees  with  respect  to 

A  1  VJ 

,  1  u 

p ,ns  D4 

Body 

Coordi nates 

Propel ler 
Coord i nates 

Real 

Imaginary 

x)04 

Th i ckness 

0.9152 

-0.2631 

9.468 

-4.01 

41 

Loading 

1.8798 

-1  .4821 

23.80 

-9.56 

35.44 

Loading  +1'* 
Thickness J 

2.7950 

-1.7452 

32.74 

-8.0 

37.02 

Bearing 

Fo  rce** 

-0.392 

-2.555 

25.705 

- 

65.30 

Net  Force*** 

-3.187 

-0.8098 

32.886 

- 

48.52 

-•'From  Hess-Smith  code  with  (+)  vertical  axis  downwards. 


■'--'•'From  DL  Propeller  Blade  Pressure  Distribution  Program  (PPEXACT)  with  (+) 
vertical  axis  upwards. 

erred  to  Propeller  coordinate  axis  with  (+)  vertical  axis  upwards. 

NOTE :  0p  is  physically  interpreted  as  being  the  blade  position  when  the 
force  on  the  body  is  a  positive  maximum  away  from  the  propeller. 

A  barogram  of  the  individual  contributions  to  the  hull  surface  force 

from  blade  thickness  and  blade  loadings  at  harmonics  of  order  0,1, 2, 3, 4  and 

5,  is  shown  in  Figure  7,  with,  in  addition,  a  comparison  of  the  net  vertical 

force  coefficient  (i.e.,  sum  of  total  surface  force  and  the  corresponding 

bearing  force)  with  the  results  of  experiments  conducted  by  the  late  Pro- 
13 

fessor  F.M.  Lewis. 

It  is  seen  that  there  is  a  difference  between  predictions  and  measure¬ 
ments  of  approx  imate  I  y  257.  In  contrast,  for  the  spheroidal  body,  the  pre¬ 
dictions  are  approximately  107  lower  than  the  measurements.  These  differ¬ 
ences  can  be  attributed  to  several  factors  but  arc  mainly  due  to  the  fact  that 
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comparison  with  Lewis'  data  is  not  strictly  valid  since  measurements  of 
the  force  include  a  combination  of  a  force  and  a  couple,  whereas  the 
result  of  these  calculations  presents  the  value  of  the  vertical  force 
only.  Therefore,  a  one-to-one  comparison  cannot,  in  principle,  be  made. 
Another  possible  source  of  uncertainty  is  the  fact  that  the  results 
have  not  yet  converged  since, when  reducing  the  mesh  size  and  increasing 
the  number  of  quadrilaterals,  the  resulting  force  keeps  changing  value. 
Furthermore,  it  is  also  possible  for  the  existing  differences  to  be 
attributed  to  a  certain  degree  to  the  fact  the  "complete"  interaction 
problem  of  the  propeller  and  hull  has  not  been  taken  into  account.  The 
"complete"  interaction  problem  involves  the  simultaneous  mutual  effect 
of  one  surface  on  the  other  carried  out  in  a  "ping-pong"  fashion  until 
stable  values  of  hydrodynamic  loadings  and  forces  on  both  surfaces  (pro¬ 
peller  and  hull)  are  established.  This  process  has  been  incorporated  in 
DL  propeller-rudder  and  propel ler-duct  programs.  The  "feedback"  of  the 
rudder  to  the  propeller  was  found  to  be  small. 

The  computing  time  in  the  present  approach  is  relatively  high  for 
practical  use,  a  fact  which  excludes  the  possibility  of  any  improvement 
of  the  scheme  by  reducing  the  mesh  size  and  increasing  the  number  of 
quadrilaterals.  A  different  procedure,  developed  under  MARAD  support, 
(mentioned  in  the  INTRODUCTION)  creates  the  hull  surface  by  a  doublet 
distribution.  This  method  is  much  more  effective  in  achieving  accuracy 
and  reduces  the  computing  time  considerably  compared  to  that  incurred 
at  present  by  using  a  hull  source  density  distribution.  The  results  of 
this  new  approach  are  also  exhibited  in  Figure  6  and  7.  The  "differences" 
for  the  spheroidal  head  have  been  reduced  to  2%,  and  for  the  Series  60 
model,  to  10%. 
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CONCLUSIONS 

A  method  based  on  rational  mechanics  is  elucidated  for  predicting 
the  hydrodynamic  forces  and  moments  acting  on  a  hull  of  arbitrary  geometry 
as  a  reaction  to  the  propeller-induced  velocity  field.  A  distribution  of 
sources  generates  the  hull  in  both  uniform  and  propeller-induced  flows,  in 
the  presence  or  absence  of  a  free  surface.  A  computer  program,  based  on 
this  procedure  and  adapted  to  the  CDC  7&00  or  Cyber  176  digital  computer, 
provides  the  strengths  of  the  sources  from  which  all  the  vertical  compo¬ 
nents  of  the  blade-frequency  prope 1 1 er- i nduced  hull  forces  can  be  found. 
Evaluation  has  been  limited  to  the  vertical  force. 

The  program  is  very  general,  i.e.,  it  is  applicable  to  any  propeller 
and  hull  of  arbitrary  form.  However,  it  requires  that  the  hull  be  sub¬ 
divided  into  a  large  number  of  quadrilaterals  in  order  to  obtain  accurate, 
reliable  results.  Since  the  propeller-induced  velocity  field  is  generally 
spatially  oscillatory  in  character,  the  dimensions  of  each  quadrilateral 
must  be  such  as  to  fit  at  least  four  of  them  into  one  cycle  of  oscillation. 
In  addition,  the  normal  velocity  to  each  quadrilateral  requires  the  evalu¬ 
ation  of  three  velocity  components  for  all  possible  combinations  of  fre¬ 
quencies  of  loading  with  those  of  the  propagation  functions.  This  part  of 
the  program  is  the  most  time-consuming  (taking  80-85%  of  the  computer 
time)  and  demands  great  attention  to  secure  convergence.  The  excessive 
number  of  quadrilaterals,  coupled  with  complications  arising  from  the 
propeller-induced  velocity,  makes  this  particular  numerical  method  expen¬ 
sive  and  rather  impractical. 

A  new  procedure  has  been  developed  under  subsequent  MARAD  support 
in  which  the  hull  is  generated  by  a  doublet  distribution.  The  doublet 
strengths  are  evaluated  by  means  of  an  integral  equation  similar  to  that 
used  in  the  present  study,  but  which  requires,  instead  of  the  three 
normal  components  of  the  propeller-induced  velocity  at  each  quadrilateral, 
the  propeller-induced  velocity  potential  (a  scalar  function)  at  the  null 
point  at  each  inner  side  of  the  panel.  In  this  way,  the  time  for  computa¬ 
tion  is  reduced  to  almost  half  that  incurred  at  present.  The  results  of 
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this  new  procedure  (exhibited  in  Figures  6  and  7)  show  much  reduced  dis¬ 
crepancies  when  compared  with  the  experimental  results.  The  "differences" 
for  the  Series  60  model  have  been  reduced  from  27%  to  107,  and  for  the 
spheroidal  body  from  11%  to  2%.  It  is  to  be  kept  in  mind,  of  course,  that 
the  force  measurements  include  a  combination  of  a  force  and  a  couple, 
whereas  the  results  of  calculations  present  values  of  vertical  force  only. 

This  method  of  doublet  distribution  is  being  employed  for  a  series 
of  surface  ship  cases  tested  by  the  late  Professor  Emeritus  F.M.  Lewis  at 
MIT  in  a  new  study,  which  will  be  reported  in  the  near  future. 
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Hu  1 1  Surface 
Water  Plane 


FIGURE  1.  PROPELLER-SERIES  60  HULL  ARRANGEMENT  WITH 
COORDINATE  AXES  WITH  RESPECT  TO  PROPELLER 
AND  TO  HULL 
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4b.  SOURCE  PANEL  DISTRIBUTION  ON  "DOUBLE"  HULL 

(SERIES  60,  CR  =  0.6,  V-SECTIONS  AFT).  ENTIRE 
HULL  PLUS  REFLECTION  IN  AN  INFINITE  FLUID, 

PORT  SIDE. 
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APPENDIX 

COMPUTER  PROGRAM 

The  program  based  on  the  foregoing  analysis  is  written  in  FORTRAN  IV 
language  and  checked  out  and  run  on  control  data  high-speed  digital  com¬ 
puters  (CDC-6600,  CDC-7600,  and  Cyber  1?6) .  The  calculations  on  CDC-6600 
and  7600  computers  were  run  with  KRONOS  and  SCOPE  operating  systems, 
respectively,  and  then  on  Cyber  176  by  making  use  of  the  NOSBE  system. 

The  final  program  of  the  "Propeller-Hull  Interaction"  problem  coded 
as  "PIHF"  is  the  result  of  a  combination  of  the  "Hess-Smi th"  program  (H-S 
code)  with  the  propeller-induced  velocity  field  program  "VFDTL"  and  the  sub¬ 
routine  designated  "FORCE"  wh i ch  implements  the  extended  "Lagally"  Theorem. 
Sketch  ft  I  exhibits  the  Flow  Diagram  connecting  the  various  subprograms  to 
generate  the  final  form  of  the  "PIHF"  code. 

Although  the  program  is  very  lengthy,  it  can  be  run  in  its  complete 
form  (i.e.,  H-S  code  +  VFDTL  +  FORCE).  It  is  advisable,  however,  to 
execute  the  program  in  stages  for  better  control  of  the  input  from  one  sub¬ 
program  to  another  and  for  keeping  track  of  the  accuracy  of  the  results  at 
each  stage. 

Initially  the  "PPEXACT"  program  is  executed  furnishing  the  propeller 
loading  information  at  all  shaft  frequencies  which  will  be  recorded  on  a 
magnet i c-tape  and  utilized  at  the  proper  time  in  the  "PIHF"  code.  Then  the 
part  of  the  "Hess-Smi th"  program  dealing  with  the  geometry  of  the  hull  form 
and  with  the  quadri laterals  into  which  the  hull  surface  is  subdivided,  will 
be  executed,  furnishing  information  about  the  direction  cosines  of  the  out¬ 
ward  directed  normal  on  each  quadrilateral,  the  coordinates  at  the  "null- 
points,"  and  the  correspond ; ng  area.  The  null  points  which  are  expressed 
with  respect  to  hull  coordinate  axes  will  be  transferred  to  the  propeller 
axes  and  saved. 

At  this  stage  the  prope 1 ler- i nduced  velocity  field  will  be  calculated 
through  the  "VFDTL"  code  at  the  null-points  previously  determined.  This 
step  is  the  most  time-consuming  one  and  should  be  done  carefully  by  plotting 
the  results  and  keeping  track  of  the  convergence. 
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These  propeller-induced  velocities  will  be  utilized  in  the  final  part 
of  the  "Hess-Smith"  program  which  deals  with  the  solution  of  the  integral 
equation.  The  solution  determines  the  strengths,  a  ,  of  the  source  distri¬ 
bution  and  these,  in  turn,  will  be  utilized  in  the  "FORCE"  coding  to  determine 
the  hydrodynamic  forces  and  moments  exerted  on  the  hull,  in  reacting  to  the 
propeller  action.  The  program,  however,  provides  meaningful  results  for  the 
vertical  component  only.  The  presence  of  the  free  surface  generates  a 
‘Momentum  flux"  through  the  waterplane  area.  The  vertical  component  of  the 
'Vnomentum  flux"  does  not  affect  the  vertical  component  of  the  hydrodynamic 
force  (see  the  Lagally  Theorem)  whereas  the  axial  and  transverse  components 
are  influenced  by  the  momentum  flux. 

Calculations  were  performed,  as  mentioned  before,  for  the  "spheroidal 
head"  and  for  a  surface  ship  of  one  of  the  Series  60,  Cn  =  0.60,  models. 

The  table  below  shows  approximate  computational  periods  running  the  program 
with  the  NOSBE  operating  system  on  the  Cyber  17&  high-speed  digital  computer; 

TABLE 

COMPUTING  TIME  IN  SYSTEM  RESOURCE  UNITS,  SRU,  ON  CYBER  176 
AND  CPU  (seconds) 

1)  Spheroidal  body  with  3-bIaded  propeller  DTNSRDC  *4118  operating  in  a 
uniform  inflow.  Body  is  subdivided  into  15*+  quadrilaterals. 


Program  CPU-Time( sec)  SRU 

PPEXACT  ( Q=0)  11.5  **1 

/H-S  (basic)  *4  15 

PIHF<VFDTL  206  738 

'H-S&Force  1 5  5*4 

238.5  (3.9*+  min.)  8U8 


2)  Series  60,  Cg=0.60,  equipped  with  5-hladed  propeller  3379*  Hull 

surface  is  subdivided  into  2C2  quadrilaterals  with  only  13*4  quadri¬ 


laterals  activated. 

PROGRAM  CPU-T i ine(  sec)*  SRU 

PPEXACT  66  236 

/H-S  (basic)  5-5  19 

PIHF  <  VFDTl(Q=0,  1,..,6)  528  1893 

^H-S  &  Force  1  5  5*' 


61*4.5  (10.2*4  min.)  2202 

If  the  CDC-7600  or  CDC-6600  computers  are  used  operating  under  KRONOS 
and  SCOPE  systems,  respectively,  then  the  times  should  be  increased 
approximately  by  a  factor  of  1.7  and  6.2,  respectively. 
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